Saturation effects in experiments on the thermal Casimir effect 
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We address three different problematic Casimir experiments in this work. The first is the classical 
Casimir force measured between two metal half spaces; here in the form of the Casimir pressure 
measurement between a gold sphere and a gold plate as performed by Decca et al. [Phys. Rev. D 
75, 077101 (2007)]; theory predicts a large negative thermal correction, absent in the high precision 
experiment. The second experiment is the measurement of the Casimir force between a metal plate 
and a laser irradiated semiconductor membrane as performed by Chen et al. [Phys. Rev. B 76, 
035338 (2007)]; the change in force with laser intensity is larger than predicted by theory. The third 
experiment is the measurement of the Casimir force between an atom and a wall in the form of the 
measurement by Obrecht et al. [Phys. Rev. Lett. 98, 063201 (2007)] of the change in oscillation 
frequency of a S7 Rb Bose-Einstein condensate trapped to a fused silica wall; the change is smaller 
than predicted by theory. We show that saturation effects can explain the discrepancies between 
theory and experiment observed in all these cases. 

PACS numbers: 42.50.Nn, 12.20.-m, 34.35. +a, 42.50.Ct 



The Casimir force is very fascinating scientifically and 
has inspired many scientists ever since Casimir published 
his classical paper [l| in 1948. It is caused by fluctuations 
in the electromagnetic fields. What is most intriguing is 
the result in the most pure geometry, the one treated by 
Casimir himself - two perfectly reflecting metal plates in 
vacuum. Here, the force is due to true vacuum fluctu- 
ations, fluctuations of the electromagnetic fields in the 
vacuum surrounding the plates. 

The interest in Casimir interactions grew very strong 
during the last decade. This increase in interest was 
spurred by the torsion pendulum experiment by Lamore- 
aux 0|, which produced results with good enough accu- 
racy for the comparison between theory and experiment 
to be feasible. This stimulated both theorists [3, \4, la, Ifl] 
and experimentalists 0, 0] and the field has grown con- 
stantly since then. Another reason for this development 
is the huge shift of general interest in the science commu- 
nity into nano-science and nano-technology where these 
forces become very important. However, the field has 
not been a complete success story. A dark cloud has 
been hovering over this field. Theory and experiment 
agree quite well for low temperatures, but at room tem- 
perature, where most experiments are performed there 
are serious deviations. Each new type of experiment 
has lead to new puzzling discrepancies between theory 
and experiment. Theorists have been forced to resort to 
phenomenological approaches to the problems, with new 
prescriptions for each new experiment. This has led to 
an unfortunate polarization of the community with those 
that are content with phenomenological descriptions on 
one side and those that want a more stringent theoretical 
treatment of the physics on the other. In this work we 
put forward what we think is the solution to the problem 
or at least the first step towards a solution. 

We address three types of experiment or experimental 



geometry: Two interacting metal plates (Gl); a semicon- 
ductor plate interacting with a metal plate (G2); an atom 
interacting with a semiconductor plate (G3). 

In all three examples the interaction energy per unit 
area, V (d), can at zero temperature be written on the 
form [l| 
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where d is the distance between the objects, k the two- 
dimensional wave vector in the plane of the plate(s), ft 
the area of a plate, and / (k, u) = is the condition for an 
electromagnetic normal mode in the particular geometry. 
The integration is along the imaginary frequency axis. 
At finite temperature the integration is replaced by a 
discrete summation over Matsubara frequencies, 
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Alternatively one may integrate along the real frequency 
axis, 
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n(o,) + l/2]ln[/(fc,u,)], (3) 



where n (lo) = [exp the distribution func- 

tion for massless bosons. This form can also be used 
at zero temperature; then the distribution function van- 
ishes. 

The force per unit area, or pressure, is obtained as the 
derivative with respect to distance, F (d) = —dV (d)/dd. 
In all three geometries there are two groups of normal 
mode, transverse magnetic (TM) and transverse elec- 
tric (TE), each with a different mode condition func- 
tion. The interaction potential is a sum of two terms, 
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Figure 1: Dispersion curves for the modes between two gold 
plates in absence of dissipation. The frequencies are in units 
of lu s , the surface plasmon frequency. The solid straight line 
is the light dispersion curve in vacuum; the dashed (dotted) 
curves are TE (TM) propagating modes; the thin solid curves 
are evanescent TM modes; the thick solid curve is the lower 
boundary for transverse bulk modes in the plates. From Ref. 



Figure 2: Casimir pressure between two gold plates. The ex- 
perimental result from Ref. [13] is shown as dots and the 
endpoints of the error bars are indicated by horizontal bars; 
the upper (lower) solid curve is the traditional theoretical 
zero (Room) temperature result; the circles are the present 
results with damping parameters 0.01, 0.1, and 1.0, respec- 
tively counting from below. 



V{d) = V™ (d) + V TE (d). In Gl and G2 the mode 
condition functions are 

TM ™> ™> 

fTE ' (k, d) = 1 - e-**W d r& E (k, u) rg E (fc, W ) , 

(4) 

where the Fresnel amplitude reflection coefficients at an 
interface between medium i and j are 

tm (i \ = gj M H (k, u) ~ £i M 7j (fc) u) 

13 1 ' 1 e ] {uj) ll {k 1 uj)+e l {u) lj {k 1 o J y W 

for TM modes (p-polarized waves) and 
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for TE modes (s-polarized waves), respectively. We have 
let the objects be of medium 1 and 2 and let the sur- 
rounding vacuum be denoted by medium 0. The gamma 
functions are 



7J (*,«) - V fc2 - £ J M ^/ c ) 5 i - 0,1,2. (7) 

Let us now study the dispersion curves for the elec- 
tromagnetic normal modes in Gl shown in Fig. Q] for 
two gold plates [lfj. This figure is valid in neglect of 
dissipation in the plate materials. The modes are propa- 
gating (evanescent) above and to the left (below and to 
the right) of the light dispersion curve. Note that there 
are no TE evanescent modes. When the system is al- 
lowed to have dissipation there are modes everywhere. 
Each original mode is replaced by a continuum of modes 



Evanescent TE modes appear and the continuum 
extends all the way down to the momentum axis. These 
modes are the cause of all the problems with the thermal 
Casimir force in this geometry. The proposed prescrip- 
tion has been to neglect the dissipation in the intraband 
part of the dielectric function but keep it in the inter- 
band part [l3|. The experimental result for the nor- 
malized Casimir pressure at 295 K is shown as dots in 
Fig. [21 The bars are the endpoints of the experimen- 
tal error bars. The upper (lower) thick solid curve is 
the theoretical result for zero temperature (295 K) cal- 
culated with Eqs. (fl]) and ([2]), respectively. We note that 
the zero temperature result agrees much better with the 
experimental result. The large negative thermal correc- 
tion comes entirely from the TE evanescent modes fl3 |. 
We will demonstrate this in more detail in a forthcom- 
ing publication [l5l |. All curves are normalized with the 
zero temperature Casimir pressure between two perfect 
metal plates, heir 2 / (240z 4 ) . We have neglected surface 
roughness effects. 

Let us now explain our view of what goes wrong in the 
theory of the thermal Casimir effect in presence of dissi- 
pation. The traditional theory relies fully on the concept 
of electromagnetic normal modes. These are assumed 
to be independent massless bosons. The possibility to 
excite one of these modes is assumed to be completely 
independent of how many modes are already excited. An 
excitation of a mode involves excitations of the charged 
particles in the system, electrons in the geometries stud- 
ied here. These are the sources of the fields. Now, the 
electrons are fermions and there is at most one electron in 
each particle state. An electron that is excited at one in- 
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Figure 3: The change in Casimir force, at 300 K, between a 
gold sphere and a silicon membrane with and without laser 
irradiation. The open squares with error bars are the experi- 
mental [T(| result. The dashed curve with open circles is the 
theoretical result without saturation effects. The solid curve 
with filled (open) circles is our present result with D equal to 
0.01 (0.1). 
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Figure 4: Fractional change in trap frequency for a Rb atom 
near a silica wall versus separation in thermal equilibrium. 
The open squares are the experimental result fl?! ]. The upper 
(lower) curve is the theoretical result including (neglecting) 
the conductivity from the few thermal carriers in the silica 
wall. The circles are our present results for the D values 
10~ 10 , 10~ n , and 10~ 12 , respectively, counted from below. 



stant of time cannot be excited again - the state is empty. 
The more modes that are excited the more difficult it is to 
excite new modes — there are saturation effects. In the 
theoretical treatment this is not taken care of. In most 
cases this fact will not cause any problems, but some- 
times it could. We think that the thermal Casimir effect 
is one such case. When dissipation is included each mode 
is replaced by a continuum of an infinite number (for an 
infinite system) of new modes. The distribution function 
diverges towards zero frequency and the saturation ef- 
fects should appear here. This is very difficult to treat 
in a strict way. We use an approximation which is very 
easy to implement. We shift the distribution function in 
Eq. Q downwards in frequency, so that it never reaches 
the point of divergence, by adding a damping parameter, 
D, 



n (oj) = [exp (H0UJ +D)-1]~ 



(8) 



The discrete frequency summation in Eq. J5]) is the re- 
sult of the poles of the distribution function that all fall 
on the imaginary axis, see Ref. Q. Our new distribu- 
tion function has its poles shifted away from the axis the 
distance D/H0 into the left half plane. The new form is 
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Each term in the summation is replaced by an integral, 
for small D values it is enough to replace only the zero 
frequency term. We will expand on this in Ref. [l5j |. The 
circles in Fig. (J2j) are the results with damping parame- 
ters 0.01, 0.1, and 1.0, respectively counting from below. 



We have used Eq. ((3} with the modified distribution 
function to get the thermal correction. 

The second geometry, G2, is a gold plate and a laser ir- 
radiated semiconductor membrane as performed by Chen 
et al. [laV They measured the change in force with 
the laser irradiation compared to without any irradia- 
tion. The results are shown in Fig. [3l The open squares 
with error bars are the experimental result. The dashed 
curve with open circles is the theoretical result for 300 
K. The deviations are clear. In this geometry it is not 
enough to neglect dissipation to get agreement with ex- 
periment. Besides, it is now the TM modes that cause 
the problems. One postulated that for the non- irradiated 
semiconductor one should completely omit the contribu- 
tion to the dielectric function from the thermally excited 
carriers. That brought the theory and experiment into 
agreement, see Fig. 10 in Ref. [la]. The solid curve with 
filled (open) circles is our saturation based result with 
D equal to 0.01 (0.1). Here we have used Eq. ^ and 
just modified the zero frequency contribution according 
to Eq. |(9|. 

The third geometry, G3, is a Rb atom near a fused 
silica wall. We study the fractional change in trap fre- 
quency versus separation when both the surroundings 
and the wall have the same temperature, 310 K. In Fig. [4] 
the experimental result []j} is shown as open squares with 
error bars. We use the same formalism as in Ref. [HI 
to find the relation between the change in trap frequency 
and the Casimir force. The upper (lower) curve is the 
theoretical result, without saturation, including (neglect- 
ing) the conductivity from the few thermal carriers in the 
silica wall. We see that also here the neglect of the contri- 
bution, to the dielectric function of the silica wall, from 
the very few thermally excited carriers brings the theoret- 
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ical result into agreement with ex per iment. This neglect 
is the postulated remedy in Ref. [la |. In this geometry, 
just as in G2, the TM modes cause the problems and it is 
not enough to neglect dissipation to get good agreement 
between theory and experiment. To include saturation 
effects we have used Eq. <[2j) and just modified the zero 
frequency contribution according to Eq. {9j. We note 
that in this experiment it is enough to have a damping pa- 
rameter as small as 10~ 10 to bring the theoretical result 
into agreement with experiment. We have assumed that 
the thermally excited carriers in the wall material have 
the conductivity 100 s" 1 (~ 10 _10 ohm _1 cm~ 1 ), which 
is the upper limit of the range given in Ref. [l8|. Us- 
ing smaller values leads to even weaker demands on the 
damping parameter. 

In summary we have proposed that saturation effects 
are responsible for the discrepancy between theory and 
experiment in several quite different Casimir geometries. 
We have treated saturation within a very simple calcula- 
tion model and demonstrated that the problems may go 
away m all cases. Other very recent theoretical models 
lgL 2(J have been proposed for the treatment of dielec- 
tric materials with a very small amount of free carriers, 
applicable to the G2 and G3 geometries. However, no 
quantitative comparison with the experiments has been 
presented. 
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